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In this paper the influence of bearing misalignment upon the dynamic characteri- 
stics of statically indeterminate rotor bearing systems is investigated. Both bearing 
loads and stability speed limits of a rotor may be changed significantly by mag n itude 
and direction of bearing misalignment. The useful theory of short journal bearings is 
introduced and simple analytical expressions, governing the misalignment problem, are 
carried out. Polar plots for the bearing load capacities and stability maps, describing 
the speed limit in terms of misalignment, are presented. These plots can be used by 
the designer to estimate deviations between calculation and experimental data due to 
misalignment effects. 


INTRODUCTION 


Bearing misalignment and initial shaft bow may considerably change the dynamic 
characteristics of a rotor bearing system. In general, bearing misalignment and shaft 
bow represent imperfections of a system with unknown amount which cannot be predicted 
exactly by the designer. Therefore, a qualitative analysis, investigating the influ- 
ence of bow and misalignment upon vibration amplitudes, dynamic bearing loads, and 
stability, etc., of a rotor, might be useful. Also significant deviations between 
experimental data and theoretical calculations of a rotor bearing system may be caused 
by bearing misalignments, for example. In this case an error estimation between ex- 
periment and calculation, considering the above mentioned imperfections, is necessary. 

It is well known from the literature that beside static and dynamic unbalance of 
a rotor a residual shaft bow may be considered as an additional external force that is 
exciting the rotor with constant amplitude independent of the speed. If a linearized 
(first order) analysis is applicable to a system then initial shaft bow and unbalance 
do not change the stability characteristics of the rotor. In other words, no interac- 
tions between self-excited shaft whirl and external excitation can be described by a 
linear approach; for example, see reference 1. Therefore, in this paper the response 
of unbalance and shaft bow is not considered, since, in a first order sense, it does 
not change the onset of instability of the rotor. 

Gasch has shown in his paper, see reference 2, that in a statically indeterminate 
rotor bearing system the static bearing loads depend upon the speed of the shaft be- 
cause of nonlinear hydrodynamic bearing characteristics. Nasuda and Hori, reference 3, 
investigated the influence of journal bearing misalignment upon the stability charac- 
teristics. They determined socalled "contour stability maps" for a two-mass and four- 
bearing rotor. 
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In this paper the theory of short journal bearings is used, see reference 4, in 
order to calculate bearing loads, bearing stiffness and damping coefficients, and sta- 
bility thresholds of a flexible rotor in terms of bearing misalignment. Both, magni- 
tude and direction of the misalignment vector of a bearing are varied and ranges of 
high and low sensitivity of the stability limit can be observed. 


NOMENCLATURE 


[A] dynamical matrix (1/s^) 

b magnitude of initial shaft bow (m) 

6 phase angle of initial shaft bow 
(rad) 

C radial bearing clearance (m) 

[Cg] bearing damping matrix (Ns/m) 

[C E ] external damping matrix (Ns/m) 

[Cj] internal damping matrix (Ns/m) 

e bearing eccentricity to clearance 
ratio 

n horizontal shaft displacement to 
clearance ratio 

{F 0 } external static load (N) 

(Fg } dynamic bearing load (N) 

{Fg 0 ) static bearing load capacity (N) 

<t> bearing load attitude angle (rad) 

[i] identity matrix 

[Kg] bearing stiffness matrix (N/m) 

[Kj] stiffness matrix due to internal 
damping (N/m) 

[Kg] shaft stiffness matrix (N/m) 

< misalignment to clearance ratio 


L bearing length (m) 

A misalignment attitude angle (direc- 

tion of misalignment with respect 
to vertical x-axis) (rad) 

m magnitude of misalignment (m) 

m , m misalignment components in x- 
x ^ and y-direction, resp. (m) 

[M] mass matrix (kg) 

U dynamic viscosity of bearing lubri- 

cant (Ns/m^) 

N + stability speed limit (Hz) 

tu angular speed of shaft (1/s) 

il> shaft center attitude angle (rad) 

R radius of the shaft at a bearing 
station (m) 


g= 

4 it C2F„ 


Sommerfeld number 


a ei genvalue ( 1 / s ) 

t time (s) 

{U} unbalance vector (kgm) 

W total weight of the rotor (N) 
{z} shaft displacement vector (m) 


5 vertical shaft displacement to 
clearance ratio 


{u} eigenvector 
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Subscripts 
B bearing 

f,l first and _last bearing station, 
respectively 

j station number (l£j£n) 

n total number of shaft stations 

o static load conditions 

x,y vertical and horizontal direction, 
respectively 


Superscripts 

T transposition of a matrix or 

column vector 

Brackets 

{ } column vector 

[ ] matrix 

( ) function of 


EQUATIONS OF MOTION 


Various definitions of bearing misalignments can be introduced for a rotor-bea- 
ring system. For example, in figure 1, lateral displacements of the bearings are mea- 
sured with respect to a line that is connecting the centers of the outer bearings 
(i.e. first and last bearing). Beside lateral bearing misalignment slope misalignment 
of the bearing axis may occur in practice but is neglected in this paper. Also slope 
deviations between shaft axis and bearing axis are not considered in this investiga- 
tion. Therefore, at least three bearings have to be present in order to get a stati- 
cally indeterminate system that exhibits misalignment effects. 

Initial shaft bow, for example, can be measured in a "non-assembled" strain-free 
configuration of the shaft where shaft centers and bearing centers at the first and 
the last bearing station coincide, see figure 1. In general, for an n-stations shaft 
the lateral misalignment vector is defined as 


{m} m {0, • . . 0,m£^ j cos j , . 

1 0, . . • 0,m£ + j sin^^ j , • 

The magnitude of misalignment m. 
tial shaft bow is given by ^ 


. m^_ | cos j ,0,. . .0| 

T 

.mi_jSinAj_j ,0, . . .0} 

for a "no-bearing" station j is defined to zero. 


( 1 ) 


Ini- 


{b(t)} = {b.cos(u)t+0. ) ,. . .0,. . .b.cos(cut+8.) , . . .0,. . .b cos(u)t+B )| 

* J J n n ( 2 ) 

|b sin(u)t+B.) ,. . .0, . . .b .sin(o>t+S. ) ,. . .0, . . .b sin(wt+8 ) } T 
» * J J n n 

in terms of the shaft rotational angle (ut, see figure 1. If internal and external dam- 
ping forces are considered and gyroscopic forces are neglected then the equations of 
motion of the rotor-bearing system can be written in the form 


[m]{£} + [c^+Oj-Hz} + [Kg+i^Hz} - 

= {F Q } + [_Kgl(b(t) } + u)^{U (t) } + {F B ({z-m},{z},u) } 


( 3 ) 
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T 

where {z} = {xj , . . .x n |y j , . . .y n } is a column vector containing lateral (vertical and 
horizontal) displacements of the shaft at n stations, see reference 1. A Taylor-ex- 
pansion of equation (3) with respect to the nonlinear bearing forces {Fg} is conver- 
gent if the static preload {F 0 } is sufficiently high compared with unbalance and 
shaft bow excitation forces. Under this assumption an expansion leads to 

[K s +K I l(z o } - {F Bo ({z o -m},{0},(J}- (F q } - {0} (4) 

and 


[M]{z'j} + [Cg+Cj+CgGm},^) ]{zj } + [Kg+KjCojJ+KgCim},^) J{Zj } - 

- [K g ]{b(t) } + ui 2 {U(t) } 


(5) 


with {z} * {z Q } + (z j (t) } being the "first order" solution of equation (3). The sta- 
tionary part {z Q } * {z o ({F o },{m}, 0 )) } of this solution is evaluated from the nonlinear 
system equations (4), with u> =* const.# 0 being the angular speed of the shaft. The 
stability of {z Q } is governed by the linear system of differential equations (5) with 
initial shaft bow and unbalance being set to zero. 

Note that first order stiffness and damping characteristics of the bearings, [Cg] 
and [Kg] respectively, depend upon the stationary solution (z 0 ({F 0 ), (m},u)) }. There- 
fore, the rotor's stability limit depends upon bearing misalignment {m} and is not 
influenced by the shaft bow {b} as long as the above expansion converges. 


BEARING CHARACTERISTICS 


Static load capacity 

The static load capacity {Fg 0 } of the bearings is evaluated by employing the 
theory of short journal bearings, see reference 4 and figure 2. For a single bearing, 
located at station j, the static bearing load components in the vertical x- and the 
horizontal y-direction can be written in the form 


C 

F„ 

Box 

< 

F 

Boy 

1 


3 

2 c 2 j 


C oj“ K j cos *j I ~Vj +K j sinA j 


r - 2 z . ^ 

oj 

2 > 2 

n .-K.sinA. | £ .-tc.cosA. 

°J 3 J j °J J J 


TT 


20-s 2 .) 3 ' 2 
^ J 


where 


03 


J 


2 2 
(5 0 j-<jCosXj) + (n o j-KjSinAj) 


( 6 ) 


(7) 


is the local bearing eccentricity ratio at station j and £ 0 j = x 9 j/Cj, n D j = y 0 j/Cj, 
and icj = mj/Cj are local static shaft displacement and bearing misalignment to bea- 
ring clearance ratios, respectively. Substituting equation (6) for each bearing into 
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equation (4) yields a nonlinear system of relations to be solved numerically for the 
static displacement vector {z Q } = {x 0 j • .x on |y 0 ] , . . .y on }^. 


Bearing stiffness and damping 
characteristics 

Once the stationary solution of the rotor at w = const.* 0 is determined for a 
given misalignment {m} , local bearing stiffness and damping matrices can easily be 
evaluated for a bearing at station j from 


4e .(1+e .) 

- 2 , 3 - 2L- 

n 2 \3 
(1 “ £ oj ) 


[K_.] = f(^p). [T.] °l _ _ ! ° J _ _ 

oj 2 -J J J o I 


and 


[c,.] [t. ] 


Bj 


3 'j L J 


-it( l+2e 2 . ) 

- . . . °J - 

2 ( 1 - e 2 .) 5/2 


ir ( 1+2e oj) 

(1-e .) 
oj' 




2e . 
oj 


f , 2 .2 

( 1-e • ) 
oj 


-4e . 
QJ 

t * 2 .2 

(1-e .) 

oj 


■2,3/2 


(1-e .) 
oj 


C V 


[T.] J 


(8) 


(9) 


CT j ] 


oj 


OJ 



5 . 

1 

-K.COSX. 1 

-n .+K.smA. 

oj 

J J 1 

1 

OJ J J 

h • 

-x.sinX. 

! E .-k.cosX. 

OJ 

J J 

1 OJ J J 


( 10 ) 


being an orthogonal transformation matrix between a local r-t-reference frame of the 
bearing j and the global yy-system, see figure 2. As an example, in figure 3a and 3b, 
respectively, dimensionless first order stiffness and damping coefficients for the 
short journal bearing are shown in terms of the Sommerfeld number. The load attitude 
angle (j> is chosen to be zero in this example, i.e. the direction of the static load 
coincides with the global (vertical) x-axis. For any other arbitrary load direction 
equations (8) and (9) have to be applied. 
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STABILITY 


After substituting the local bearing stiffness and damping matrices into the glo- 
bal matrices of equation (5) the corresponding eigenvalue problem 

[A - <jl] {u} = {0} (11) 

can be solved for a as a function of the rotor angular speed w and the bearing misa- 
lignment vector {m}. The dynamical matrix [A] is of dimension 2n x 2n and is given by 


[A({m},cu) ] 


[0] I 

1 

[I] 

1 

-Im]” 1 [Kg + K^U))]- j 

-Cm] _1 [c e + Cj.]- 

-[M] -1 [K B ({m},u)] j 

-[M] _1 [C B ({m>,u))] 


( 12 ) 


EXAMPLE 


Figure 4 shows a symmetric five stations rotor with three short journal bearings 
and two major masses mounted on a flexible shaft. The data for this rotor are listed 
in Table 1 . The center bearing of the system may be misaligned with respect to the 
outer bearings by and my in a vertical and horizontal direction, respectively. The 
numerical calculations for this example were carried out on the main frame computer 
CDC-CYBER 170/730 of the Technical University of Vienna, Austria. The nonlinear equa- 
tion solver ZSPOW and the eigenvalue solver EIGRF of the IMSL-computer library were 
used. 

Figure 5 shows the static load ratio of the center bearing in terms of the shaft 
speed. The horizontal misalignment my is kept zero while the vertical misalignment m x 
is varied. It can be seen that each curve in figure 5 approaches the stability thres- 
hold at a certain speed of the rotor. The bearing load gradient is high at low speeds 
and decreases with increasing speed. The stability threshold is increasing with in- 
creasing misalignment above 12 um and below -5 ym in this example. In general, the 
center bearing load ratio in figure 5 is increasing with incrasing speed except with- 
in a small range of misalignments between (-30 ym £ £ -15 ym) at low speeds. A si- 

milar diagram can be drawn for the outer bearing load ratios. Note that these bearing 
load characteristics might change with changes in the shaft bending stiffness and the 
bearing clearances. Also, in practice, the misalignment attitude angle \ is unknown 
and different from zero. Therefore, the influence of both mx # 0 and my * 0 upon bea- 
ring loads and stability limits has to be studied. 

Figure 6 shows polar plots for the bearing load ratios with the rotor speed kept 
constant at N=200 Hz. Curves of constant magnitudes of misalignment are drawn in a 
plane where the polar angle indicates the direction of misalignment (attitude angle \ ) 
and the radius is the bearing load ratio, i.e. at constant radii the bearing loads are 
constant and at constant polar angles the ratio my/mx is constant. At zero misalign- 
ment the load ratios of both center and outer bearings are constant. Hence, the curves 
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m=0 represent circles with constant bearing loads. For increasing magnitudes of misa- 
lignment the bearing load ratios change substantially with the direction of misalign- 
ment X. For example, if a misalignment of 10 pm is present in the system, (i.e. 40% 
of the center bearing clearance or 1/24000 of the outer bearing span) then the center 
bearing load ratio varies between 0.2 at X=25° and 1.09 at X=225°. A corresponding 
change for the outer bearing can be seen on the right hand side of figure 6. 

It is to expect that for misalignment directions that cause low bearing loads in 
the system the stability speed limit of the rotor is decreased also. This can be ob- 
served from figure 7, showing a polar plot of the stability limits of the present ro- 
tor. The polar angle in this diagram again indicates the direction of misalignment X 
while the radius now is the stability speed limit. Curves of constant misalignment 
magnitudes are drawn in the figure. For low misalignments up to m=12 pm and for atti- 
tude angles of about X=20° and X=220° the stability limits tip down significantly. For 
example, if m=10 pm, then the stability limit of the rotor changes between N + =275 Hz 
and N + =365 Hz. Within the ranges of (180° £ X £ 250°) and (0£m£12 pm) a high sensiti- 
vity of the stability speed limit with respect to the direction of misalignment can 
be observed. 

Stability maps and polar plots as discussed above give a complete information to 
the designer and make it possible for him to estimate misalignment effects upon the 
stability thresholds and the bearing load capacities of a given rotor. 


CONCLUSION 


The theory of short journal bearings makes it possible to formulate the misalign- 
ment problem of statically indeterminate rotor hearing systems in a simple and easy 
manner. The numerical effort to carry out bearing loads and stability limits in terms 
of misalignment is comparatively low even for a multi--station rotor supported by a 
multi-bearing system. There are a number of good and reliable algorithms available 
in various computer libraries to solve the governing equations of the problem. Since 
magnitude and direction of bearing misalignment are never exactly known in a real sy- 
stem the above presented theory is sufficiently accurate for practical applications 
and can be used by the designer in order to make a good estimation on whether or not 
a certain amount of bearing misalignment might dynamically endanger or damage a rota- 
ting machinery. Also the above results may explain deviations between experimental 
data and predictions where misalignment effects were not included. 
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DIMENSIONLESS BEARING STIFFNESS AND DAMPING CHARACTERISTICS 
FOR THE SHORT JOURNAL BEARING 
(Load attitude angle $ » 0°) 


STIFFNESS 



Figure 3a. 



SOMMERFELD NUMBER 
Figure 3b. 


SYMMETRIC ROTOR BEARING SYSTEM 
WITH MISALIGNMENT 



Figure 4. 


CENTER BEARING LOAD RATIO VS. 
SHAFT SPEED FOR VARIOUS VERTICAL 
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Figure 5. 
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BEARING LOAD RATIOS Fg/W IN TERMS OF MISALIGNMENT ATTITUDE 
ANGLE A AND MISALIGNMENT MAGNITUDE M 
(Rotor speed N * 200 Hz) 



Figure 6. 


STABILITY SPEED LIMITS N IN 

TERMS OF MISALIGNMENT ATTITUDE DATA OF THE ROTOR BEARING SYSTEM 

ANGLE A AND MISALIGNMENT MAGNI- AS SHOWN IN FIGURE 4 

TUDE M 



Total mass of the rotor 
Shaft bending stiffness 
Outer bearing span 


16.8 kg 
3350 Nm 2 
240 am 


BEARING DATA 


CENTER 


OUTER 


Length 

Diameter 

Radial clearance 
Lubricant vicosity 


25 mm 
25 um 


12.5 mm 
15 um 


0.01 Ns/m 2 


Table 1. 


Figure 7. 



